Anderson impurity in a correlated conduction band 
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We investigate the physics of a magnetic impurity with spin 1/2 in a correlated metallic host. 
Describing the band by a Hubbard Hamiltonian, the problem is analyzed using dynamical mean- 
field-theory in combination with Wilson's nonperturbative numerical renormalization group. We 
present results for the single-particle density of states and the dynamical spin susceptibility at 
zero temperature. New spectral features (side peaks) are found which should be observable ex- 
perimentally. In addition, we find a general enhancement of the Kondo scale due to correlations. 
Nevertheless, in the metallic phase, the Kondo scale always vanishes exponentially in the limit of 
small hybridization. 

PACS numbers: 71.27.+a, 75.20.Hr 
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The Anderson model M has been successfully applied 
in the past to describe the physics of a magnetic impurity 
embedded in a conducting host. Extensive theoretical 
studies of this particular many-body problem led to con- 
siderable insight as well as progress in the development of 
new methods 0. The most thoroughly analyzed case is 
that of an impurity in a noninteracting conduction band 
with a constant density of states. Several properties have 
been established by Wilson's numerical renormalization 
group ||,[| and by the Bethe ansatz ||. Most impor- 
tantly, a new many-body energy scale Tk (the Kondo 
temperature) arises, which is exponentially small in the 
limit of vanishing hybridization. In addition, it has been 
demonstrated || that below this temperature the system 
can always be understood as a "local Fermi liquid" with 
strongly renormalized quasiparticles. The single-particle 
spectrum was shown to exhibit a generic three peak struc- 
ture consisting of two atomic levels and a quasiparticle 
resonance of width Tk- Based on this model, a number of 
experimental results for dilute impurities in metals have 
been explained successfully, including measurements of 
the resistivity, the magnetic susceptibility and the spe- 
cific heat. For a review see ref. @]. 

It has become clear, however, that the single impu- 
rity Anderson model is somewhat too simplified and that 
qualitatively different types of physical behaviour are 
possible when a more general Hamiltonian is considered. 
One very important characteristic of real materials is the 
interaction among the conduction electrons. This aspect 
is usually neglected, mostly for technical reasons, i.e. to 
simplify the investigation. If taken into account, we ex- 
pect, at least, a renormalization of the model parame- 
ters. Our work will focus on the question whether, in 
addition, qualitatively new physics is possible. An ex- 
perimental realization frequently cited in this context is 
the cuprate system Nd2_ x 0e x CuO4 ||], a concentrated 
impurity system, where the energy scale of low tempera- 
ture heavy fermion behaviour is apparently incompatible 
with the standard Kondo picture. 



Models with a single impurity embedded in a cor- 
related host were studied already within several ap- 
proaches. Perturbative calculations in a slave boson rep- 
resentation by Khaliullin and Fulde || yielded a renor- 
malization of the effective Kondo coupling. Very simi- 
lar results were obtained by Tornow at al. ]ic[ | within a 
non-crossing approximation. Furthermore, in the limit of 
high dimensions and using a variational treatment, Davi- 
dovich and Zevin |ll[] found a qualitative change of the 
behaviour of the Kondo temperature Tk- According to 
these authors, above some intermediate value of the con- 
duction band interaction, Tk is no longer exponentially 
small at vanishing exchange coupling. In our work, we 
will discuss this issue in detail. Finally, in the case of one 
dimension, Phillips and Sandler ]T3] and also Schiller and 
Ingersent jl3| ] represented the interacting host as a Lut- 
tinger liquid, which makes a renormalization group treat- 
ment possible. Among other results, they found that in 
some region of parameter space an unquenched local mo- 
ment may survive down to the lowest temperatures. 





FIG. 1. Left: Anderson impurity (square) coupled to one 
site of a correlated lattice (circles). The hybridization is taken 
to be purely local. Right: effective two-impurity model with 
noninteracting band orbitals (empty circles) of energy and 
hybridization Vi. 

These studies already indicate the competition of sev- 
eral effects: (i) The conduction band correlations may 
change the density of states (DOS) of the conduction 
band, (ii) A repulsive on-site interaction will reduce the 
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hybridization of the impurity level, (iii) The conduc- 
tion electrons will become increasingly polarized, thus 
enhancing the effective spin coupling of the impurity mo- 
ment. In the following, using Wilson's nonperturbative 
numerical renormalization group, we will analyze which 
one of these factors dominates. 

Our Hamiltonian consists of a spin- 1/2 impurity em- 
bedded in an interacting host (see also fig. |l|) represented 
by a one-band Hubbard model |14|: 



H = - ]T (t^ - e c 5 l3 ) c| CT Cja + U B J2 n *T n H 

ija i 
a 

Note that the impurity hybridizes with a single conduc- 
tion band orbital, which in the following will be denoted 
as the i = Hubbard site. 

We will be interested in the case of half filling, which 

- assuming a bipartite lattice and next-neighbour hop- 
ping (tij = t) only - is equivalent to e c = —Ub/2 and 
ef = —(7/2. The calculations in our paper are restricted 
to the paramagnetic phase of the host. A controlled ap- 
proach |l5|,|l6| to correlated lattice problems is possible 
in the limit of large coordination number Z — > oo, scaling 
the hopping matrix elements as t = -K^ . In our treatment 
we will use the Bethe lattice and take the noninteracting 
half-bandwidth D = 2t* = 1 as the unit of energy. It 
should be emphasized, however, that the choice of the 
lattice is merely motivated by calculational convenience 
and should have no qualitative effect on the results JlTj . 

It was pointed out in Q that integrating out all 
the band fermions except those on the i — Hub- 
bard site (the so-called "cavity method") yields an ef- 
fective action which has the same retarded part iS c ff,ret = 

— JJ dr dr' cl a Qq 1 (t — t')co ct (t') as the pure Hubbard 
model. As a result, the system can be described as a 
two- impurity model (see fig. [l]) with an effective noninter- 
acting bath defined by a hybridization function A c (cj) — 
ir J2 P | V^, | 2 c5(c^ — e p ) . In our calculation we will therefore 
follow a two-step procedure: First, we solve the Hubbard 
model in dimension d — > oo using Wilson's numerical 
renormalization group |^],|| as in G^j. In the paramag- 
netic phase considered here this leads to a Mott transition 
at a critical interaction strength U B * 2.92. Correlations 
strongly influence the structure of the DOS; close to the 
transition point an effective narrow band Jf8| is formed 
by the quasiparticle resonance. 

In the next step, we add the /-impurity. The combined 
system is then again treated using NRG, this time with- 
out the self-consistency loop (the modification of the ef- 
fective bath due to the impurity is 0(1/ N) and can there- 
fore be neglected in the thermodynamic limit). The band 
correlations enter via the previously determined DOS and 
the c-site interaction Ub- In our calculations we have 



confined ourselves to the metallic regime Ub < 2.92 and 
zero temperature. 

First, we present results for the single particle spec- 
tra pff c )(co) = — i ImGy( c )(w) of the impurity / and the 
Hubbard site i = 0. Considering pf (fig. 2a), we ob- 
tain a three-peak structure in the spectrum already for 
7 = and small hybridization A = 355- • This may be 
attributed to a narrowing of the effective band, leading to 
resonances at finite energy ]l8| ]. Upon increasing U, these 
peaks are reduced and weight is shifted to the atomic lev- 
els which for large interaction can be found at u> ~ ±7/2. 
In addition, the height of the quasiparticle peak is sig- 
nificantly reduced. Luttinger's theorem |lj| which states 
that p f (0) is pinned at its U = value is therefore found 
to be not valid in the case of an interacting conduction 
band. 

Furthermore, we notice that the width of the quasipar- 
ticle resonance is almost independent of U, in contrast 
to the situation with Ub — 0. This already indicates 
a strong enhancement of the Kondo scale due to band 
correlations, as will be discussed in more detail below. 




FIG. 2. Spectral densities pf (a) and p c (b) at zero temper- 
ature, A = 0.1 and Ub = 2.6. For comparison, we also show 
the c spectral density at vanishing hybridization A, equivalent 
to the DOS of the Hubbard model. 

The corresponding results for p c are shown in fig. 2b. 
For any finite A a hybridization gap is formed. This 
is already the case at U = and Ub = (not shown 
here) and persists for finite interactions, indicating that 
the system is a Fermi liquid (a nonvanishing self-energy 
at lu = would smear out the gap). The Fermi liquid 
picture is independently supported by the fact that the 
fixed point of the NRG and its leading irrelevant eigen- 
operators are unchanged compared to the noninteracting 
case; see also [poi. 
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A quantity which is more easily accessible experimen- 
tally is the (longitudinal) dynamic susceptibility, defined 
as the response of the impurity spin to a local magnetic 
field 



(2) 



X(u;)=i / dte iut ([SI(t),Sf(0)]). 



Within the NRG formalism it is convenient j2l| to cal- 
culate the imaginary part x"^) directly and to obtain 
the real part y i a Kramers-Kronig transformation. 

In particular, the static susceptibility is then given by 
Xo = x'(0). In the following we will focus on the spin 
relaxation function 



S(u,) = 



(3) 



In a first step, we consider its behavior with increasing 
band correlations for fixed A, U (see fig. 3a). 



U B =0.0 
U B =2.0 
U B =2.5 
U B =2.8 
U B =2.9 




FIG. 3. Local spin relaxation function for fixed U = (a) 
and Ub = 2.9 (b). In both cases, the hybridization has the 
value A = 0.01. Note that due to the normalization a value 
of 1 at ui = would be consistent with the Shiba relation. 

In this plot we employed a normalization of x"(w) sug- 
gested by the Shiba relation [J22|]23 



(4) 



derived for Ub — p3[ ]. For the noninteracting band 
this relation is indeed satisfied (with an error of less than 
10% due to the NRG procedure). With increasing Ub 
systematic deviations arise, indicating that the general 
proof p3] based on Ward identities breaks down for an 



interacting conduction band. The lineshape of x"(<x>) also 
depends on Ub ■ For a weakly correlated band we obtain 
a single elastic peak, while close to the metal-insulator 
transition at U B = 2.92 two additional inelastic side 
peaks arise. They indicate that in this case the Kondo 
singlet is formed at an energy scale which lies outside the 
effective band. The width of the remaining elastic peak 
is determined by the effective bandwidth of the Hubbard 
model. 

An increase of the impurity interaction U (fig. 3b), 
leads to a suppression of the elastic peak and to a shift 
of the inelastic peaks (corresponding to a slight reduction 
of the singlet binding energy) . 

For the real part x'(w), some typical results are shown 
in fig. [|: Already at a weak band interaction U b (when 
the Hubbard DOS is well approximated by the nonin- 
teracting one) the static susceptibility xo 1S strongly re- 
duced. 




FIG. 4. Real part of the dynamic susceptibility at A = 0.1 
and U = 1.0. 

Of particular interest is the hybridization dependence 
of the low energy scale, the Kondo temperature Tk , at 
intermediate to strong band interaction. While there 
is agreement on the fact that a small Ub enhances the 
effective Kondo coupling j^-[ll],^o| but still leads to an 
exponentially vanishing Tk at small hybridization A, it 
was found in ref. ||ll]| that above an intermediate U b the 
Kondo temperature varies linearly in A. We will now 
consider this issue in detail. 

We define Tjc to be equal to the binding energy of the 
local singlet, which is given by the position of the maxi- 
mum in x"(w). Note that this definition also applies to 
the case of a noninteracting impurity (U=0). With in- 
creasing A we observe a crossover from an exponential 
to a power law behavior Tk ~ A. The crossover point 
depends on Ub and is proportional to the effective band- 
width D c tf. For very small A, the Kondo temperature 
always varies as InTjf ~ —U/A. In contrast to ref. |Tl| ] 
we therefore find an exponentially small Tk at any Ub, as 
long as the host is metallic. The discrepancy may be due 
to the approximate variational method used in ref. pd| ]. 

At a finite band interaction, Ub can lead to a non- 
monotonic behavior of Tk] see fig. |. The increase at 
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small Ub can be attributed to the local interaction on 
the site i = while the decrease close to the MIT is 
due to band narrowing. As Ub — > £4iit, the Kondo 
scale approaches a finite limiting value, indicating that 
even in the paramagnetic insulator the local impurity is 
screened! We can understand this by considering the 
effective hybridization "seen" by the /-impurity 
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Lu + iO+ -A c (uj+) 
In the insulating host A c = and therefore 
A /H ~ V 2 5((j). 



(•5) 



(6) 



In this case, the impurity couples exclusively to the i = 
site, the singlet is purely local and no Kondo many- 
particle physics is possible. 
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FIG. 5. Kondo temperature as a function of the band in- 
teraction. 

In conclusion, we analyzed a model describing an An- 
derson impurity in a correlated band. The limit of large 
coordination numbers made a treatment within the dy- 
namical mean-field-theory possible. Using the nonper- 
turbativc numerical rcnormalization group, we solved the 
corresponding effective two-impurity model and obtained 
the one particle spectra and the dynamic susceptibility. 
We found that the system is always a Fermi liquid as 
long as the host is metallic. Band correlations lead to a 
strongly enhanced Kondo scale, indicating that the dom- 
inant effect of Ub is to increase the spin polarization of 
the conduction electrons. Nevertheless, Tk remains ex- 
ponentially small as a function of hybridization. This is 
consistent with a Fermi liquid picture of the Hubbard 
host where the Kondo screening of the impurity is due 
to fermionic quasiparticles instead of bare electrons. In 
the spectral quantities, a change of the lincshape and the 
formation of side peaks is observed close to the Mott tran- 
sition. This is explained by a narrowing of the effective 
conduction band. 

Of course it would be desirable to compare our results 
with experiments on systems that can actually be consid- 
ered as dilute. One possibility might be to perform ESR 



measurements on rare earth systems, where impurities 
can be introduced into a correlated host in a controlled 
way. Here the main experimental signal (the absorption 
X"{oj)) could be directly related to our model calcula- 
tions. 

In future calculations we will study the effects of finite 
temperature, different fillings and antiferromagnetic or- 
der on our findings. We will also extend our analysis to 
thermodynamic and transport properties. 

The authors would like to thank T.A. Costi, H.- 
A. Krug von Nidda, A. Schiller and S. Tornow for useful 
discussions. 
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